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Abstract 

We analyze the direct detection signals of a toy model consisting of a Dirac dark matter 
particle which couples to one Standard Model fermion via a scalar mediator. For all scenar¬ 
ios, the dark matter particle scatters off nucleons via one loop-induced electromagnetic and 
electroweak moments, as well as via the one-loop exchange of a Higgs boson. Besides, and 
depending on the details of the model, the scattering can also be mediated at tree level via the 
exchange of the scalar mediator or at one loop via gluon-gluon interactions. We show that, 
for thermally produced dark matter particles, the current limits from the LUX experiment 
on these scenarios are remarkably strong, even for dark matter coupling only to leptons. We 
also discuss future prospects for XENONIT and DARWIN and we argue that multi-ton xenon 
detectors will be able to probe practically the whole parameter space of the model consistent 
with thermal production and perturbativity. We also discuss briefly the implications of our 
results for the dark matter interpretation of the Galactic GeV excess. 


1 Introduction 

Dark matter direct detection experiments have witnessed in the last years an impressive improvement 
in sensitivity, reflected in particular by the recent breaking of the zeptobarn-scale barrier in the 
scattering cross section by the LUX experiment [T] . Moreover, the prospects for future improvements 
are bright; in the near future the XENONIT experiment [2] will extend the reach in cross section 
by about one order of magnitude and, in the longer term, a multi-ton xenon detector, such as 
DARWIN |3], might even reach the yoctobarn scale. 

The excellent sensitivity of experiments has allowed to rule out already some well motivated 
models, notably those where the dark matter is in the form of Weakly Interacting Massive Particles 
(WIMPs) which couple at tree level to the Standard Model via the exchange of the Z-boson, such as 
the sneutrino in the Minimal Supersymmetric Standard Model |1]. On the other hand, some other 
scenarios have suppressed tree-level couplings to the nucleons and can naturally evade detection at 
current experiments. This is the case when the dark matter interaction with the nucleus is mediated 
by the Higgs boson, and hence suppressed by the small Higgs coupling to the proton constituents 
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(for example, when the dark matter is a singlet scalar [5l E], an inert SU{2)l donblet [71 El |9] a 
hidden vector [10] or a singlet antisymmetric tensor m)- Alternatively, this is also the case when 
the dark matter only conples to the leptons or to the heavy qnarks, where the tree-level conplings to 
the nucleons simply vanish (as occurs, for instance, in the radiative seesaw model [12]) • Nevertheless, 
even if the dark matter-nucleon tree-level coupling vanishes, such a coupling is necessarily generated 
via quantum effects, due to the non-vanishing gauge and Higgs coupling interactions of all the 
Standard Model fermions. Given the excellent sensitivity of current instruments, such quantum 
effects could induce observable signals at experiments. The strength of the signal is, however, very 
model dependent and should be investigated case by case. 

In this paper we focus on scenarios with a Dirac singlet fermion as dark matter candidate which 
couples to a Standard Model fermion via a scalar mediator. The model contains only three free 
parameters, namely the dark matter mass, the mediator mass and the Yukawa coupling strength, 
which can be hxed in terms of the other two parameters from requiring that the dark matter popula¬ 
tion in our Universe is generated via thermal freeze-out. In this paper we a present a comprehensive 
analysis of the scattering process of dark matter particles with nucleons for couplings to any Stan¬ 
dard Model fermion. The scattering arises at tree level in the case of the hrst generation quarks and 
radiatively for other fermions, namely via the one-loop exchange of a Higgs, a Z-boson and a photon 
for any fermion [niEllISlESllITlEHlliniEniEI], or via a box-diagram with two external gluons 
for colored fermions [22l [23l [23| . We then calculate the scattering rate with the nucleons induced by 
this effective interaction, which we hnally confront with the upper limit reported by LUX, as well 
as with the projected reach of the XENONIT and DARWIN experiments. 

The paper is organized as follows. In section we describe the model. Then, in section we 
discuss the effective dark matter-nucleon interaction induced in the model, for dark matter coupling 
to any Standard Model fermion. In section [^ we calculate the corresponding scattering rate with 
nucleons for each of these possibilities and we confront the predicted rate to the present limits from 
the LUX experiment and to the projected sensitivities of future experiments. Finally, in section]^ 
we present our conclusions and in Appendix [^ our method to calculate the event rate in LUX. 

2 Singlet Dirac dark matter with a charged mediator 

We consider a simple extension of the Standard Model by a Dirac fermion y, singlet under the 
Standard Model gauge group and which constitutes our dark matter candidate, and a charged scalar 
particle rj, which mediates the interaction between the dark matter and a Standard Model fermion 
/. Furthermore, in order to guarantee the stability of the dark matter particle, we impose a discrete 
Z 2 symmetry under which y and r] are odd, while all the Standard Model particles are even. The 
most general from for the Lagrangian reads: 

n _ n \ r \ r I Y»fermion i /»scalar /-i \ 

- >^SM + + Mnt +Mnt * UJ 

Here, £sm is the Standard Model Lagrangian, which includes a potential for the Higgs doublet d), 

V = -I- |Ai(<hi<I))^, while and £ri are the parts of the Lagrangian involving just the dark 

matter particle or the scalar mediator and which are given by: 

= X0X - m^XX and (2) 

, (3) 
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where the covariant derivative depends on the gange qnantnm nnmbers of rj, which will be 
specihed below. 

On the other hand, describes the Ynkawa interaction between the dark matter particle, 

the scalar mediator and the Standard Model fermion /, which can be either a right-handed fermion 
singlet fa G {mr, Cj?}; * = l!2,3 being a generation index, or a left-handed fermion donblet 
Jl G Then, the interaction Lagrangian reads 

4:r”" = -!/')'x/R + h.c.. (4) 

for conpling to a right-handed fermion and 

^fermion f-r (-yVoX^^L - VV-XSl + h-c. for conpling to , and 

= —yvxjL + ii-c. = < , + (5) 

\ -y vIxul - y rilxdi + h.c. for conpling to Ql 

for conpling to a left-handed fermion. Notice that specifying the fermion representation hxes the 
gange qnantnm nnmbers of the scalar mediator. 

Lastly describes the interaction between the scalar mediator and the Standard Model 

Higgs boson and reads: 


_y 3 (<|)t<j))(j^t^) for conpling to/ij, and 

_y 3 (<|)t<j))(j^t^) _ for conpling to fi ■ 


( 6 ) 


These interaction terms only affect the phenomenology in a signihcant way when the conplings are 
0{1). Therefore, in what follows we will neglect these interactions and we will simply set these 
conplings to zero. This choice implies in particnlar that, for the scenario where the dark matter 
conples to left-handed fermions, the two components of the scalar donblet rj have a common mass, 
which we denote by rrir^. 

Provided that the conpling y is large enongh, the dark matter candidate y is kept in thermal 
eqnilibrinm in the Early Universe with the plasma of Standard Model particles. For large parts of the 
parameter space, the annihilation process which is most relevant for the eqnilibration is XX //, 
where / is the Standard Model fermion that conples to the dark matter particle via the interaction 
given in Eq. Q or (|^. The thermally averaged annihilation cross section at the time of dark matter 
deconpling is, for Dirac fermions, well approximated by the velocity independent part of the cross 
section, which reads: 


(av) 


~ av\ 



(7) 


In this expression, Nc = 3 for conpling to qnarks and Nc = 1 for conpling to leptons. Moreover, 
in scenarios where the dark matter particle conples to a fermion donblet fi, the total annihilation 
cross section is given by the snm of the annihilation cross sections into the two states of the fermion 
donblet. As is well known, for Dirac dark matter the total annihilation cross section for the process 
XX ff s-f the time of freeze-ont shonld take the valne (crn)thermai — 4.4 x 10“^®cm^/s in order 
to match the dark matter density in onr Universe, = 0.1198 ± 0.0015, as measnred by the 

Planck satellite pg. 

In this work we also consider the possibility that the scalar mediator rj is close in mass to the dark 
matter particle x, fo which case the calcnlation of the dark matter abnndance mnst be modified dne to 
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the presence of coannihilation processes. Namely, if rriri ^ 1.2 rriy^, the temperature at the time of dark 
matter decoupling can be large enough to keep r] also in thermal equilibrium with the plasma, which 
efficiently depletes the number of dark matter particles through the processes XV qg, —)■ gg, 
and gug\ W~^g, if x couples to quarks, and through xV e“ 7 , gg"' —>■ 77 , 7 Z, ZZ, hh+hh“, or 
7oh- iy’'' 7 , if X couples to leptons. In order to fully take into account all coannihilation channels 
in the calculation of the relic density, we use micrOMEGAs [26], interfaced with FeynRules |2^ and 
CalcHEP [2H| for the numerical solution of the corresponding Boltzmann equation. 

In our numerical analysis we will focus in the well motivated scenario where the Dirac dark 
matter particle is produced via thermal freeze-out in the early Universe. This condition Exes one 
of the parameters of the model in terms of the remaining two, e.g. the Yukawa coupling y = 
j/thermai(?^x’Then, the parameter space of the model is spanned by only two parameters, for 
example, the dark matter mass and the relative difference between the dark matter mass and 
the mediator mass, (m,, — 

3 Effective Lagrangian for dark matter-nucleon scattering 

The interactions of a Dirac dark matter particle with a nucleon can be described by the Feynman 
diagrams shown, up to one-loop in perturbation theory, in Fig. where we include for conciseness 
only one representative diagram of each class. Panel (a) corresponds to the tree-level interaction 
mediated by the exchange of the scalar g. Panels (b) and (c) show penguin diagrams mediated by 
the photon (or the Z-boson) and by the Higgs boson, respectively, which arise at the one loop level 
in all models of Dirac dark matter with a charged scalar mediator, regardless of the choice of the 
Standard Model fermion. Finally, panel (d) shows a box-diagram with two external gluons, which 
arises also at the one loop level when the dark matter couples to a colored scalar mediator. Let us 
discuss separately the effective Lagrangian that arises in each case. 

If the dark matter couples at tree level to a first generation quark, the dominant contribution 
to the scattering cross section is the tree level process XQ XQ with the f-channel exchange of the 
scalar mediator g, depicted in Fig. panel (a). Analogously, the dark matter antiparticle x scatters 
off quarks, with the scalar being in this case exchanged in the s-channel. After a Fierz rearrangement 
of the corresponding matrix elements, these diagrams give rise to a vector interaction of x with the 
quark, which in turn translates into an effective vector interaction of the dark matter particle with 
the nucleon of the form: 


CeS,tree = XYX ( 8 ) 

and which leads, in the non-relativistic limit, to a spin-independent interaction. Here, the effective 


couplings for protons and neutrons, fy^ 

and read 



(2\v 

for coupling to Ur 

and = 1 

f Ay 

for coupling to Ur 

f(p) — } \ 

J ytree S 

for coupling to dR 

1 2 Ay 

for coupling to dR 

Is Ay 

for coupling to {urYl, 

) 1 

Is Ay 

for coupling to (ml, di) 


with = y^ j [8 (m^ vector-current conservation, the spin-independent tree-level 

contribution is exactly zero for coupling to second- and third-generation quarks, and is of course 
also zero for dark matter coupling solely to a lepton. The same diagrams also give rise to an axial 
vector interaction of the form f^^eeXYYxYtiX^Q^ which leads to a spin-dependent interaction. The 
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Figure 1: Exemplary Feynman diagrams contributing to the dark matter-nucleon scattering: (a) tree-level 
scattering, (b) one-loop exchange of a photon or Z boson, (c) one-loop exchange of a Higgs boson, and (d) 
scattering off gluons at one loop. 


coupling constants of the axial vector and the vector interaction turn out to be comparable in size in 
this model, therefore, given that the experimental sensitivity to the spin-independent interaction is 
orders of magnitude better than to the spin-dependent one, this contribution can be in most instances 
safely neglected. Nevertheless, it will also be included in our numerical analysis for completeness. 

On the other hand, and regardless of the choice of the Standard Model fermion, the Yukawa 
interactions Eq. Q and Eq. (|^ necessarily induce at the quantum level the scattering of dark 
matter particles with nucleons, due to the penguin diagrams depicted in panels (b) and (c) of Eig. 
The first diagram induces electromagnetic moments for the dark matter particle, y, the most relevant 
ones for Dirac dark matter being the magnetic dipole moment /i^ and the charge radius b^, defined 
by the effective dark matter-photon Lagrangianj^ 

. ( 10 ) 


The electromagnetic moments jj,^ and are obtained by matching the coefficients of this effective 
Lagrangian to the results of an explicit calculation of the relevant loop diagrams. For dark matter 
coupling to a right-handed Standard Model fermion and taking the limit where the transferred 
momentum in the dark matter-nucleus interaction goes to zero, the result reads0 




-QfeNcy'^ 
327r^m^ 
-QfeN^y‘^ 


2 / A ^/2 

(8A^ -|- (9/i + 7e — 5) A — 4e(3/i -|- e — 1)) arctanh 


( 11 ) 


fi € — 1 


+ (8fi - 8e + 1) log + 4 ^4 + 1 


( 12 ) 

Here, Qf is the charge of fn in units of e (such that e.g. Qf = —1 for coupling to cr), Nc = 3 (1) 


Mhis diagram also leads in addition to an anapole moment, which corresponds to the effective operator 
This operator is suppressed in the non-relativistic limit by the square of the dark matter velocity 
and is always subdominant in this model. This contribution, however, is also included in our numerical analysis for 
completeness. 

^We have used FeynCalc [H] for parts of the computation. 
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for coupling to quarks (leptons) and 


A = /r + (e — 1) — 2/i(e + 1) , 


(13) 


with /i = m^/m^ and e = (see also [IE]). For dark matter coupling to a left-handed Standard 

Model fermion doublet, /i^ and are simply given by the sum of the corresponding expressions for 
both components in the doublet. 

The expression for presents, when nif rurj, a infrared divergence oc log , which 


enhances the scattering rate for light fermions. Note, however, that Eq. (12) was derived in the limit 


where the transferred momentum goes to zero and is therefore not valid in dark matter scenarios 
where the fermion mass constitutes the smallest energy scale in the scattering process. More specif¬ 
ically, when the dark matter scatters off xenon nuclei, the typical transferred momentum is ~ 50 


MeV, therefore for couplings to electrons and first generation quarks Eq. (12) is not valid. For these 


cases, and in order to cut-off the divergence, we replace in our numerical analysis the fermion mass 
by 50 MeV. On the other hand, and in contrast to the charge radius, the magnetic dipole moment 
has a hnite value for m/ -C 

Analogously, a dark matter interaction with the Z boson arises at the one-loop level, leading to 
an effective vector interaction of the form: 


^eS,Z = fv^z Xl^X , 


(14) 


where 

f(p) _ Gpaz An) _ Gpaz 

Jv,z — ^ ) Jv,z — ^ 

Here, az is an effective form factor which, for dark matter coupling to a right-handed fermion, is 
given by 


AIr) 


Qj '7 — 


TiN, 


c2/^e 


Ibvr^ 


^ (- 

2 V/i 


+ 


1 A 


AV2 


arctanh 


AV2 


/i -f- e 


(16) 


with A dehned as in Eq. (13). For the case of dark matter coupling to a left-handed fermion doublet 
one has —and Notice that the dark matter effective coupling 

to the Z boson scales as making this contribution subdominant in comparison to the 

photon exchange, except for the interesting scenario where the dark matter couples to tn or (ti, bi)- 
The penguin diagram shown in panel (c) of Fig. [^induces a coupling of y to the Standard Model 
Higgs boson h, which in turn induces the effective dark matter-nucleon interaction 


Higgs /s, Higgs XX -VA , 


(17) 


where 


AN) 

7s,Higgs 


—V^GFm^mNfN^ 3e 

ml 327r2A^/2 

2 (e^ -|- /X (/i — 1) — e (1 -|- 2/i)) arctanh 


AV2 


jj. € — 1 


+ A^/^ 1^2 + (/i 


e) log 





(18) 
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with the Higgs-nucleon coupling ~ 0.345 [U] j As the coupling of the fermion / to the Higgs is 
proportional to rrif, this contribution is clearly subdominant for dark matter coupling to any of the 
light Standard Model fermions. In the case of dark matter coupling to a third-generation quark, we 
find that the contribution from the Z exchange is always much larger, rendering the Higgs exchange 
diagram subdominant for all possible cases, although it is also included in our numerical results. 
An additional contribution to the effective dark matter-Higgs interaction arises when the quartic 
couplings A 3 and A 4 in Eq. @ are non-vanishing. However, we find numerically that, even for 
quartic couplings 0{1), the Higgs interaction gives a subdominant contribution to the dark matter 
scattering off nucleons. 

Finally, if the dark matter couples to quarks, it can scatter off gluons via one-loop box diagrams, 
such as the one shown in panel (d) of Fig. This generates an effective scalar interaction of y with 
nucleons, which reads 




efF, gluon 


= /SL„ XXiViV, 


(19) 


The effective coupling has been calculated for the case of Majorana dark matter interacting 

with a colored mediator in [221 ESI ES] • By explicitly writing down the amplitudes for all relevant 
diagrams, it is straightforward to check that /g^iuon exactly the same for Majorana and for Dirac 
dark matter. Hence, we do not repeat here the rather cumbersome expression for the effective scalar 
coupling, and we refer instead to these works. Also, following the detailed discussion in m , we 
include the gluon contribution to the effective dark matter-nucleon Lagrangian only for the case 
where dark matter couples to a heavy quark, i.e. c, b or t. 


4 Event rate for dark matter-nucleus scattering and exper¬ 
imental limits 


The effective Lagrangians given in Eqs. 14, 17, 19) describing the interaction of dark matter 


with nucleons, as well as the effective dark matter-photon interaction defined by Eq. (10), induce 
the elastic scattering of dark matter off nuclei. Taking into account all these contributions, the 
differential scattering cross section of dark matter off a nucleus with mass ttla and spin Ja which 
recoils with energy Er, reads, including only the leading terms in an expansion in 1/Er and 1/u^: 


dcr 

d^ 




2^2 


1 


rriA 




[Fsi {Er)Y 


I^Al^lmA Ja + 1 


:x 

TTU 

rriA 


' 3Ja 

(/'■"T l^si (£«)] 


[^dipole (^h)I* 


( 20 ) 


where 


is the effective scalar dark matter-nucleus coupling, defined as 


= Z 


Ap) 

•I S, gluon 


, Ap) 

' J ytree 


, Ap) 
^ Jv,z 


-eb^- 


2m^ 


'j+iA-Z) (/^ 


(n) r{n) An) 

gluon ' •'Vjtree ' J V,Z 


( 21 ) 


^The effective dark matter interactions to the Higgs and to the Z-boson have also been presented in |3D] in the 
limit ^ 0 for a model where the dark matter couples to tR. Our results in that limit agree with the expressions 
given in that work. 
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For Fsi{Eji), we use the Helm form factor given in |3T], while for the form factor Fdipoie (-F_r) 
associated to the dipole-dipole scattering we follow [32] . 

Finally, we translate the scattering cross section into an event rate in a xenon experiment and 
we compare the predicted rate with the current limit by the LUX experiment pQ, which provides 
the strongest sensitivity for the model discussed in this paper. It is important to remark that one 
can not simply map the limits on cxgj as they are published by the corresponding collaborations onto 
limits on the parameter space of the model considered in this work, due to basically two reasons: 
hrst, the scattering diagrams involving the exchange of a photon only induce a dark matter coupling 
to protons, but not to neutrons, leading to a (potentially large) source of isospin violation in the 
scattering. Secondly, due to the long-range force induced by the magnetic dipole moment of dark 
matter, the differential scattering cross section given in Eq. (20) has a different dependence on the 
recoil energy Er compared to the case of standard spin-independent or spin-dependent interactions. 
In our numerical analysis, we derive the corresponding upper limits arising from the null-result of 
LUX, fully taking into account the isospin violation as well as the non-standard dependence of the 
differential scattering cross section with Er. We provide the details of our method of calculating 
the event rate in LUX in Appendix 

We will also calculate the regions of the parameter space which can be probed by the future 
experiments XENONIT [2], already under construction since 2013, and DARWIN, a design study 
for a next generation multi-ton dark matter detector [3]. Following [33|, we assume that XENONIT 
will have a reach in a factor of 40 better than the present LUX limit for dark matter masses above 
~ 100 GeV. For the model discussed in this work, this translates into a rescaling of the upper limits 
on the Yukawa coupling y arising from the LUX data by a factor of 40^^^ ~ 2.5. Here, we refrain 
from deriving prospects for XENONIT for small dark matter masses, i.e. ruy. < 30 GeV, as these 
will depend on so far unknown details of the detector performance near threshold. For DARWIN, 
we consider the most optimistic scenario where the reach in cross section reaches ~ 10“^® cm^ |3], 
thereby approaching the irreducible background of coherent neutrino-nucleus scattering, and we 
derive prospects by rescaling our limits on the Yukawa coupling obtained from the LUX data by a 
factor of 1800^/'^ ~ 6.5. As for XENONIT, we do not consider prospects for low-mass dark matter. 
For simplicity we also assume that DARWIN will consist only of xenon, although currently both 
liquid argon and xenon are discussed as possible target materials. 

Before presenting the actual limits on the model, let us hrst discuss the relative importance 
of the different contributions to the scattering rate, which were introduced in Section and were 
schematically shown in Fig. [T] Which of these processes gives the largest contribution to the event 
rate in a direct detection experiment depends on which Standard Model fermion the dark matter 
particle couples to, as well as on the dark matter and mediator masses. In Figs, [^andj^we show, 
separately for each coupling scheme, which processes dominate the scattering in LUX, for given 
dark matter and mediator masses, and m^. The left panels of Fig. correspond to a scenario 
where the dark matter particle couples to right-handed up-quarks (top left), right-handed down- 
quarks (central left), or to the hrst generation of left-handed quarks (bottom left). As expected, 
in these cases the tree-level scattering of dark matter off the nucleons dominate the event rate, 
rendering all one-loop processes irrelevant. On the other hand, as discussed in Section the tree- 
level contribution is negligible if dark matter couples to quarks of the second or third generation, 
or to leptons. In all these scenarios, the one-loop processes dominate the event rate in LUX (and 
other direct detection experiments), as shown in Figs. and More precisely, the scattering via 
photon exchange, which leads to the charge radius and magnetic dipole moment operators, is the 
dominant contribution to the cross section for large parts of the parameter space of each scenario. In 
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0.1 


tree-level scattering 
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1 0.1 
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Figure 2: Dominant contributions to the dark matter scattering rate in LUX, for coupling to right-handed 
up-type quarks (upper panels), right-handed down-type quarks (middle panels), and left-handed quark 
doublets (lower panels). 


particular, the charge radius operator typically leads to the largest event rate for < 1 TeV, while 
the magnetic dipole moment dominates for > 1 TeV, although the concrete values depend on 
the coupling scheme and on the mass splitting between the mediator r] and the dark matter particle 
X- This behavior can be qualitatively understood from the dimensionality of each effective operator 
describing the coupling of dark matter to photons (see Eq. (10)), and which is dimension hve for 
the case of magnetic dipole interactions, while it is dimension six for the charge radius. Hence, the 
corresponding coefficients of these operators scale differently with the dark matter mass, namely as 
l^x ~ ^ respectively. 

On the other hand, if the dark matter particle couples to the top-quark, the scattering via the 
exchange of a. Z boson becomes relevant, and generically dominates over the photon exchange for 
masses below ~ 5 — 10 TeV, as shown in the right upper and right bottom panels of Fig. As 
mentioned in Section this is due to the dependence of the effective coupling az with the mass of 
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Coupling to or (z/g, e^) 


Coupling to or 


Coupling to Tjf or 




Figure 3: Dominant contributions to the dark matter scattering rate in LUX, for coupling to leptons. 


the fermion coupling to from Eq. (16) it follows that for m/ this effective coupling scales 

as az oc {mf/m^Y, making it relevant only for the case of x coupling to tji or ipL^tL). Finally, we 
hnd that the dark matter-nucleon scattering induced by the one-loop coupling of dark matter to 
gluons is only relevant for rather small dark matter masses and small mass splittings rriri/m^, for 
dark matter coupling to a heavy quark, as illustrated in the corresponding panels in Fig. 

We hnally determine the regions of the parameter space which are excluded by the LUX ex¬ 
periment for each of these scenarios, as well as the regions that can be probed with the future 
XENONIT and DARWIN experiments, hxing the value of the Yukawa coupling at each point to 
the value i/thermal required to reproduce the observed dark matter abundance via thermal freeze-out. 
The results for the case of dark matter coupling to quarks are shown in Fig. and for coupling 
to leptons in Fig. The dotted black contour lines show the value of the Yukawa coupling which 
leads, at that point of the parameter space, to the observed dark matter abundance — 0.12. 

Besides, for each of the possible scenarios, the region of the parameter space corresponding to low 
mass and small mass splitting is theoretically inaccessible since the coannihilations suppress the relic 
abundance of x below the observed value. This region is shown in Figs. and as dark grey. On 
the other hand, large masses and/or large mass splittings correspond to Yukawa couplings y > -\/47r, 
where our perturbative calculation loses validity]^ It is worthwhile noticing that this requirement 
sets an upper limit on the dark matter mass which approximately reads uidm ^ 20 TeV for coupling 


to quarks and judm ^ 15 TeV for coupling to leptons. 

The left panels in Fig. correspond to scenarios where x couples to hrst generation quarks, 
namely to right-handed up-quarks (top panel), right-handed down-quarks (central panel), and to 
the hrst generation left-handed quark doublet (lower panel). As apparent from the plots, practically 
the whole parameter space for a thermal relic compatible with our perturbativity condition is already 
excluded by the LUX experiment, with XENONIT being able to constrain dark matter masses up to 
~ 10 TeV, thereby probing the whole parameter space of the model. These strong direct detection 
bounds result from the unsuppressed tree-level coupling of the dark matter particle x to light quarks, 
cf. Eq. (g. 

On the other hand, for coupling to quarks of the second generation, the tree-level coupling 
is negligible and the limits dominantly arise from the dark matter electromagnetic moments, cf. 
Fig-i Consequently, the constraints from LUX are considerably weaker, as apparent from the 


^This condition is motivated by the requirement that the expansion parameter ay = for higher-order 

corrections should be smaller than one. 
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Figure 4: Regions of the parameter space for Dirac dark matter particles with scalar mediator excluded 
by the LUX experiment, assuming coupling to right-handed up-quarks (upper panels), right-handed down- 
quarks (middle panels), and left-handed quark doublets (lower panels), and assuming that the dark matter 
density of the Universe is generated by thermal freeze-out (the necessary value of the Yukawa coupling 
is indicated by the dotted black contour lines). The dark gray shaded regions in the lower left corners 
of each figure are theoretically inaccessible, due to efficient coannihilations which deplete the dark matter 
density below the observed value, while the light gray shaded regions violate our perturbativity condition 

ythermal Y \/dvT. 


second column of Fig. |^for dark matter coupling to cr (top panel), sr (central panel), and (c^, sl) 
(lower panel). Nevertheless, it is worth pointing out that despite the suppression of the one-loop 
diagrams, LUX already excludes dark matter masses up to ~ 100 — 500 GeV, depending on the mass 
splitting of X and r] as well as on the precise coupling scheme of dark matter. Most interestingly, in 
these scenarios there are exciting prospects for the upcoming direct detection experiments to probe 
large parts of the viable parameter space. For example, a DARWIN-like experiment would be able 
to cover practically the whole parameter space of thermally produced dark matter, compatible with 
our perturbativity condition, despite the fact that all direct detection processes occur only at the 
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one-loop level. Furthermore it is worthwhile remarking that direct detection experiments have an 
enhanced sensitivity in the case of mass-degenerate scenarios, i.e. rrir^/m^ ~ 1, as apparent from the 
hgures. This larger sensitivity follows from the parametric dependence of the charge radius operator 

oc — 1), which leads to an enhancement of the scattering cross section as the mass 

spectrum becomes more and more degenerate. 

Lastly, for dark matter coupling to third-generation quarks, the constraints on the parameter 
space are presented in the right panels of Fig. In the case of coupling to tn, shown in the 
upper right panel, we restrict our analysis to > rrit — 175 GeV, which ensures that the process 
XX tt is kinematically allowed during freeze-out in all of the considered parameter space. In 
principle, a thermal relic can also exist for rriy. < rrit, however in that case the relic density would 
instead be set by one-loop annihilations or two-to-three processes, changing the phenomenology 
of the model quite dramatically. For an extensive discussion of the model in that part of the 
parameter space we refer to [30]. As argued in section if the dark matter particle couples to 
either tji or {tL,bi), the scattering rate in direct detection experiments is typically dominated by 
the Z exchange contribution, while it is dominated by the electromagnetic moments for coupling to 
the right-handed bottom quark. Then, for coupling to the right-handed top quark, LUX excludes 
masses up to 200 — 400 GeV, depending on the mass splitting of rj and y, while the constraints 
for coupling to fo/j (central panel) are somewhat weaker due to the absence of a large effective dark 
matter coupling to the Z boson. Again, it is important to point out that XENONIT and eventually 
DARWIN will probe a significantly larger region of the parameter space than LUX, and will be 
able to rule out, or hopefully discover, dark matter particles with masses up to several TeV. Similar 
conclusions also hold for the case of y coupling to the third generation quark-doublet, shown in the 
lower right panel of Fig. 

Direct search experiments are subject to astrophysical and nuclear uncertainties which might 
modify the conclusions presented above. Hence, and in order to robustly exclude a region of the 
parameter space or to conhrm a possible signal, it is important to exploit the complementarity 
with indirect and collider searches, which rely on very different physical processes, and which are 
affected by different systematic uncertainties. A thorough investigation of the complementarity 
among the different search strategies for Dirac dark matter is beyond the scope of this paper and 
here we will limit ourselves to briefly comment on the qualitative features of the indirect and collider 
searches. If the colored mediator is sufficiently light, it might be directly produced at the LHG via 
strong interactions. If, moreover, the mass difference between the scalar mediator and the dark 
matter particle is sizable, the subsequent decay produces a characteristic signal consisting in two 
jets plus missing transverse momentum, which has been searched for in the framework of simplified 
supersymmetric models jSH [35l |36l |3^. In contrast, when the mass difference is small, the jets 
might be too soft to be detected and a different search strategy should be pursued, for instance 
the search for monophoton [3H] or monojet events [39] plus missing transverse momentum. Finally, 
for very heavy mediators, only the dark matter particle can be pair-produced at the collider. In 
this regime, dark matter signals could be detected in searches for monojet, monophoton or mono- 
W/Z boson plus missing transverse momentum |10lHTlll2]. Present collider limits constrain some 
regions of the parameter space with < 200 GeV, although the precise reach of the limits crucially 
depends on the details of the model. Additional constraints from colliders follow from electroweak 
precision observables, notably the Z-boson decay width |13], and which exclude the existence of 
scalar mediators with mass < 40 GeV. 

On the other hand, the strongest indirect detection limits on this class of scenarios stem from the 
non-observation of an excess of gamma-rays correlated to the direction of dwarf galaxies. A recent 
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analysis by the Fermi-LAT collaboration from a stacked analysis of 15 dwarf spheroidal galaxies 
exclndes ttidm ^ 100 GeV for dark matter conpling to the 6 -qnark H, which is competitive, and 
for some choices of parameters better, than the LUX limits shown in Fig. Additional constraints 
on the model can be derived from gamma-ray observations of the Milky Way center, althongh this 
search is affected by large, and still not fnlly understood, backgrounds. Interestingly, some authors 
have claimed evidence for a gamma-ray excess in this region ng 061 HU EHl sni ED], which can 
be interpreted in terms of dark matter annihilations. Dirac dark matter with a charged mediator 
constitutes an excellent particle physics framework to explain the excess (in contrast to Majorana or 
scalar dark matter, which have annihilation cross sections at the Galactic center which are helicity- 
and velocity-suppressed). The mass ranges favored by the dark matter interpretation of the galactic 
center excess (GGE) for the annihilation final states selected in [51], and which point to an annihi¬ 
lation cross section of the order of the thermal value, are shown in Fig. as vertical hatched bands. 
It follows from the hgure that direct detection experiment can test, for some channels, this interpre¬ 
tation of the excess, especially for couplings to hrst generation quarks, to sr or to cr. One should 
bear in mind that the astrophysical uncertainties in the extraction of the dark matter parameters 
from the GeV excess are still large. Therefore, the value of the annihilation cross section can be 
significantly lower than the thermal value and accordingly the impact of the LUX limits. 

Finally, we turn to scenarios in which the dark matter particle only couples to one of the right- 
or left-handed leptons of the Standard Model. As discussed in Section in this scenario, the tree- 
level interactions vanish, however the one-loop diagrams induce significant scattering rates in direct 
detection experiments. In fact, many of the qualitative features of the scenario where y couples 
to heavy quarks also apply to this case. The parameter space for the various Dirac dark matter 
scenarios with coupling to leptons are shown in Fig. where again the region of the parameter 
space incompatible with thermal production and perturbativity are shown, respectively, as dark 
and light gray shaded regions. For coupling to leptons, the region of parameter space for which no 
thermal relic exists due to strong coannihilations is smaller compared to the case of dark matter 
coupling to quarks, since for leptophilic dark matter the relevant coannihilation channels (such as 
rjr]^ — 77 ) are suppressed by powers of the electromagnetic coupling, while for coupling to quarks, 
the corresponding processes are induced by the strong coupling between the scalar rj and the gluons. 
Besides, in the case of dark matter coupling to leptons, for ~ ~ 45 GeV the coannihilation 

process 77 I —>■ // is resonantly enhanced due to the diagram where the Z boson is exchanged in the 
s-channel. Hence, around that dark matter mass, the region in parameter space where no thermal 
relic can exist due to coannihilations is extended to higher mass splittings. 

The region of parameter space excluded by LUX, as well as the range of masses which will be 
accessible to XENONIT and DARWIN are very similar for the six different scenarios shown in Fig. 
For all mass splittings between the mediator 7 and the dark matter particle y, LUX excludes dark 
matter masses between 8 GeV < < 100 GeV, and XENONIT will be sensitive up to ~ 500 

GeV. Moreover, as noted already before, the charge radius operator 6 ^ is enhanced for small mass 
splittings, leading to an increased sensitivity of the direct detection experiments for —>■ 1 . 

Moreover, from these figures it follows that a multi-ton detection experiment like DARWIN will 
practically cover the whole parameter space of these scenarios, again despite the fact that all direct 
detection signals are one-loop suppressed. 

Scenarios with Dirac dark matter coupling to leptons can also be probed in collider or indirect 
search experiments, which give complementary limits on the parameter space of the model. The 
charged scalar mediator can be pair-produced via the Drell-Yan process in colliders and subsequently 
decay producing a characteristic signature of opposite-sign, same-flavor leptons, plus missing trans- 
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Figure 5: Same as Fig. but for dark matter coupling to right-handed charged leptons (upper panels) 
and to left-handed lepton-doublets (lower panels). 


verse momentum. This search has been undertaken in the context of slepton searches at the LHC 
by the ATLAS [36] and CMS collaborations [37], and at LEP II by the ALEPH, DELPHI, L3 and 
OPAL collaborations ISI. and exclude some regions of the parameter space with ^ 100 GeV (see 
[53] for a compilation of the limits of the parameter space of a Majorana dark matter particle that 
couples to a lepton and a scalar mediator, and which can also be applied to Dirac dark matter, since 
the production mechanism of scalar mediators and their decay modes are identical in both cases). 
Finally, the precise measurement of the Z-boson decay width at LEP [13] excludes mediators with 
mass < 40 GeV. 

Stringent limits on this scenario also follow from indirect searches. The annihilation into leptons 
produces an additional component in the cosmic positron flux which is severely constrained by the 
AMS-02 data on the positron fraction [SU |35], dehned as the flux of positrons divided by the total 
flux of electrons plus positrons, as well as the data on the positron flux [56]. More specihcally, the 
non-observation of the sharp feature in the positron fraction produced by the hard leptons produced 
in these annihilation channels translates into strong limits on the parameters of the model. From 
the limits on the cross section derived in Ea one can exclude < 120 GeV for coupling to cr, 

< 70 GeV for coupling to fiR and < 30 GeV for coupling to tr. The more conservative limits 

from [SH] exclude rriy. < 100 GeV and < 40 GeV for couplings to cr and ^r, respectively. Strong 

constraints on these channels can also be derived from dwarf galaxy observations, which exclude 

< 100 GeV [44] . 
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5 Conclusions 


In some scenarios, the dark matter particle only couples to the nucleons via quantum effects, with 
a scattering rate which is expected to be very small. On the other hand, the excellent sensitivity of 
current direct search experiments might allow to identify such rare events over the environmental 
and instrumental backgrounds, thus opening the possibility of probing a class of models which were 
not accessible to previous experiments, such as those where the dark matter couples only to leptons 
or to heavy quarks. 

In this paper, we have investigated this exciting possibility in the framework of a toy model 
where the dark matter is a Dirac fermion, singlet under the Standard Model gauge group, which 
couples to one of the Standard Model fermions via a Yukawa coupling with a scalar mediator. For 
reasonable values of the parameters of the model, the observed dark matter abundance ~ 0.12 

can be explained via thermal freeze-out. In this model, the scattering rate with nucleons receives 
various contributions, which depend on the choice of the Standard Model fermion. For couplings 
to first generation quarks, the dark matter interacts with the nucleons at tree level resulting in 
a comparatively large scattering rate. On the other hand, for couplings to any Standard Model 
fermion, the dark matter interacts with the nucleon via one-loop penguin diagrams mediated by the 
photon, the Z boson or the Higgs boson. Lastly, for couplings to any colored fermion, the dark 
matter interacts with the nucleon also via one-loop box diagrams with two external gluons. 

Remarkably, current experiments have reached the necessary sensitivity to probe these quantum- 
induced dark matter interactions with nucleons. More specifically, assuming thermal production for 
the dark matter population in our Universe, the current LUX data probe the parameter space of all 
scenarios considered in this paper, posing in some cases very stringent constraints on the parameter 
space of the model. For instance, the LUX data exclude dark matter masses below ~ 3 TeV for 
coupling to first generation quarks, masses between ~ 10 — 100 GeV for coupling to the right- 
handed charm or strange quarks, and masses also between ~ 10 — 100 GeV for coupling to leptons 
(of any generation and chirality). The future XENONIT experiment will continue closing in on the 
parameter space of the model and will be able to probe dark matter masses up to several hundreds of 
GeV, while in the longer term, a multi-ton xenon detector such as DARWIN will be able to exclude 
masses as large as a few TeV. Direct search experiments will then have a high sensitivity to this 
class of scenarios, even for multi-TeV dark matter masses, which are difficult to probe at the LHG. 
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A Derivation of the limits from the LUX data 

We recall from Section]^ that the limits on the spin-independent scattering cross sections published 
by the LUX collaboration cannot be applied to our scenario, due to the large isospin violation and, 
especially, due to the different dependence of the differential cross section with the recoil energy, as 
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a result of the scatterings mediated by a photon. In this appendix we briefly describe our method to 
translate the LUX null results into limits on the Dirac dark matter model considered in this paper. 


Starting from the differential dark matter-nucleus scattering cross section, Eq. (20), the number 


of expected recoils per kilogram target mass and per keV recoil energy can be straightforwardly 
calculated, the result being: 


dRi 


iiP® 


dE 


R ruN^rUx Jv^uEr) 


dn vfE{v) 


dcr 


dE 


R 


with Vrain{Eji) = 


ImNiEji 

^hred 


( 22 ) 


In this expression, the index i refers to a given xenon isotope with mass mAr- and mass fraction 
Pq ~ 0.3 GeV/cm^ is the local dark matter density (of which p©/2 is made up by x and Po/2 by 
x), and fE{v) is the one-dimensional speed distribution of dark matter particles in the Earth’s rest 
frame. The latter is derived from an isotropic Maxwell-Boltzmann distribution in the Galactic rest 
frame with velocity dispersion uq = 220 km/s, using a mean Earth’s velocity of ve = 233.3 km/s. 

The total expected scattering rate in the LUX reads [59] 
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Here ^{Ee) is the expected number of photoelectrons (PE) for a nuclear recoil with energy E^ which 
we extract from Fig. 4 of [T] by taking the intersections of the red solid line with the contours of 
constant Er. The actual number of photoelectrons n is then given by a Poisson distribution with 
mean u{Ee)] following the LUX collaboration, we conservatively do not model any recoils below 
= 3 keV. We take into account the single-photoelectron resolution of the photomultipliers by 
assuming that the observed Si signal is given by a Gauss distribution with mean n and standard 
deviation y/napMT, for which we take ctpmt = 0.37 [50]. The acceptance e (Si) of the LUX detector 


is given by the red dashed line in Fig. 1 of [T], and the additional factor 1/2 in Eq. (23) takes into 
account that we dehne only the region below the mean nuclear recoil band (red solid line in Fig. 4 
of m) as the signal region for dark matter recoils; see [61] for an in-depth discussion of this choice. 
Finally, S'™™ = 2 PE and = 30 PE are the minimal and maximal number of photoelectrons 
considered in the LUX analysis. 

In the signal region, LUX observed one event after running with an exposure of 85.3 x 118.3 kg ■ 
days. Making no assumptions about the number of expected background events, we obtain a 90% 
G.L. upper limit of iVmax = 3.89 expected recoil events, which then can be used to set limits 
on the normalization of the dark matter scattering cross section. As a cross-check we compare 
our limits on an isospin-conserving spin-independent interaction with those published by the LUX 
collaboration [1]: depending on the dark matter mass, our limits are more conservative by a factor 
1.5-2.5. 


References 

[1] LUX Collaboration, D. Akerib et. ai, First results from the LUX dark matter experiment 
at the Sanford Underground Research Facility, Phys.Rev.Lett. 112 (2014), no. 9 091303, 
|arXiv: 1310.8214|. 

[2] XENONIT collaboration, E. Aprile, The XENONIT Dark Matter Search Experiment, 
Springer Proc.Phys. C12-02-22 (2013) 93-96, |arXiv: 1206.6288|. 


16 










[3] DARWIN Consortium, L. Baudis, DARWIN: dark matter WIMP search with noble liquids, 
J.Phys.Conf.Ser. 375 (2012) 012028, |arXiv: 1201.2402], 

[4] T. Falk, K. A. Olive, and M. Srednicki, Heavy sneutrinos as dark matter, Phys.Lett. B339 
(1994) 248-251, |hep-ph/9409270|. 

[5] J. McDonald, Gauge singlet scalars as cold dark matter, Phys.Rev. D50 (1994) 3637-3649, 

|hep-ph/0702143j. 

[6] J. M. Cline, K. Kainnlainen, P. Scott, and C. Weniger, Update on scalar singlet dark matter, 
Phys.Rev. D88 (2013) 055025, |arXiv: 1306.4710|. 

[71 N. G. Deslipande and E. Ma, Pattern of Symmetry Breaking with Two Higgs Doublets, 
Phys.Rev. D18 (1978) 2574. 

[8] L. Lopez Honorez, E. Nezri, J. F. Oliver, and M. H. Tytgat, The Inert Doublet Model: An 
Archetype for Dark Matter, JCAP 0702 (2007) 028, |hep-ph/0612275|. 

[9] D. Majnmdar and A. Ghosal, Dark Matter candidate in a Heavy Higgs Model - Direct 
Detection Rates, Mod.Phys.Lett. A23 (2008) 2011-2022, |hep-ph/0607067|. 

[10] T. Hambye, Hidden vector dark matter, JHEP 0901 (2009) 028, |arXiv:0811.0172|. 

[11] O. Gata and A. Ibarra, Dark Matter Stability without New Symmetries, Phys.Rev. D90 
(2014), no. 6 063509, jarXiv: 1404.04321. 

[12] E. Ma, Verifiable radiative seesaw mechanism of neutrino mass and dark matter, Phys.Rev. 
D73 (2006) 077301, |hep-ph/0601225|. 

[13] P. Agrawal, S. Blanchet, Z. Ghacko, and G. Kilic, Flavored Dark Matter, and Its Implications 
for Direct Detection and Colliders, Phys.Rev. D86 (2012) 055002, jarXiv: 1109.3516|. 

[14] Y. Bai and J. Berger, Fermion Portal Dark Matter, JHFP 1311 (2013) 171, 

|arXiv: 1308.0612|. 

[15] Y. Bai and J. Berger, Lepton Portal Dark Matter, JHFP 1408 (2014) 153, |arXiv: 1402.66961. 

[16] S. Ghang, R. Edezhath, J. Hntchinson, and M. Lnty, Leptophilic Fffective WIMPs, Phys.Rev. 
D90 (2014), no. 1 015011, |arXiv: 1402.7358]. 

[17] P. Agrawal, Z. Ghacko, and G. B. Verhaaren, Leptophilic Dark Matter and the Anomalous 
Magnetic Moment of the Muon, JHFP 1408 (2014) 147, |arXiv: 1402.7369|. 

[18] P. Agrawal, B. Bated, D. Hooper, and T. Lin, Flavored Dark Matter and the Galactic Center 
Gamma-Ray Excess, Phys.Rev. D90 (2014), no. 6 063512, |arXiv: 1404.1373|. 

[19] J. Kopp, L. Michaels, and J. Smirnov, Loopy Constraints on Leptophilic Dark Matter and 
Internal Bremsstrahlung, JCAP 1404 (2014) 022, |arXiv: 1401.6457|. 

[20] Z.-H. Yu, X.-J. Bi, Q.-S. Yan, and P.-F. Yin, Tau Portal Dark Matter models at the LHC, 
arXiv:1410.3347, 


17 



[21] A. DiFranzo, K. I. Nagao, A. Rajaraman, and T. M. Tail, Simplified Models for Dark Matter 
Interacting with Quarks, JHEP 1311 (2013) 014, |arXiv: 1308.2679|. 

[22] M. Drees and M. Nojiri, Neutralino - nucleon scattering revisited, Phys.Rev. D48 (1993) 
3483-3501, |hep-ph/9307208|. 

[23] J. Hisano, K. Ishiwata, and N. Nagata, Gluon contribution to the dark matter direct detection, 
Phys.Rev. D82 (2010) 115007, |arXiv: 1007.2601]. 

[24] P. Gondolo and S. Scopel, On the sbottom resonance in dark matter scattering, JCAP 1310 
(2013) 032, larXiv: 1307.4481]. 

[25] Planck Collaboration, Planck 2015 results. XIII. Cosmological parameters, 

arXiv:1502.01589, 

[26] G. Belanger, F. Bondjema, A. Pnkliov, and A. Semenov, micrOMEGAs 3: A program for 
calculating dark matter observables, Comput.Phys.Commun. 185 (2014) 960-985, 

|arXiv: 1305.0237|. 

[27] A. Allonl, N. D. Gliristensen, G. Degrande, G. Duhr, and B. Fnks, EeynRules 2.0 - A complete 
toolbox for tree-level phenomenology, Comput.Phys.Commun. 185 (2014) 2250-2300, 

|arXiv: 1310.1921|. 

[28] A. Belyaev, N. D. Gliristensen, and A. Pnkhov, CalcHEP 3.4 for collider physics within and 
beyond the Standard Model, Comput.Phys.Commun. 184 (2013) 1729-1769, 

|arXiv: 1207.6082|. 

[29] R. Mertig, M. Bohm, and A. Denner, PEYN CALC: Computer algebraic calculation of 
Eeynman amplitudes, Comput.Phys.Commun. 64 (1991) 345-359. 

[30] A. Knmar and S. Tnlin, Top-flavored dark matter and the forward-backward asymmetry, 
Phys.Rev. D87 (2013), no. 9 095006, |arXiv: 1303.03321. 

[31] J. Lewin and P. Smith, Review of mathematics, numerical factors, and corrections for dark 
matter experiments based on elastic nuclear recoil, Astropart.Phys. 6 (1996) 87-112. 

[32] T. Banks, J.-F. Fortin, and S. Thomas, Direct Detection of Dark Matter Electromagnetic 
Dipole Moments, arXiv: 1007.5515. 

[33] P. Gnshman, G. Galbiati, D. McKinsey, H. Robertson, T. Tait, et. ai. Working Group Report: 
WIMP Dark Matter Direct Detection, arXiv: 1310.8327, 

[34] CMS Collaboration, S. Ghatrchyan et. ai. Search for new physics in the multijet and 
missing transverse momentum final state in proton-proton collisions at \/s= 8 TeV, JHEP 
1406 (2014) 055, | arXiv: 1402.47701. 

[35] ATLAS Collaboration, G. Aad et. ai. Search for sguarks and gluinos with the ATLAS 
detector in final states with jets and missing transverse momentum using y/s = 8 TeV 
proton-proton collision data, JHEP 1409 (2014) 176, |arXiv: 1405.7875|. 


18 






[36] ATLAS, G. Aad et. ai, Search for direct production of charginos, neutralinos and sleptons in 
final states with two leptons and missing transverse momentum in pp collisions at \/s = 8 
TeV with the ATLAS detector, JHEP 1405 (2014) 071, |arXiv: 1403.5294|. 

[37] CMS, V. Khachatryan et. ai, Searches for electroweak production of charginos, neutralinos, 
and sleptons decaying to leptons and W, Z, and Higgs bosons in pp collisions at 8 TeV, 
Eur.Phys.J. C74 (2014), no. 9 3036, |arXiv: 1405.7570]. 

[38] ATLAS Collaboration, G. Aad et. ai. Search for new phenomena in events with a photon 
and missing transverse momentum in pp collisions at ^/s = '8 TeV with the ATLAS detector, 

arXiv:1411.1559, 

[39] ATLAS Collaboration, G. Aad et. ai. Search for pair-produced third-generation sguarks 
decaying via charm guarks or in compressed supersymmetric scenarios in pp collisions at 
^/s = 8 TeV with the ATLAS detector, Phys.Rev. D90 (2014) 052008, |arXiv: 1407.0608|. 

[40] ATLAS Collaboration, G. Aad et. ai. Search for new phenomena in final states with an 
energetic jet and large missing transverse momentum in pp collisions at y/s = 8 TeF with the 
ATLAS detector, arXiv: 1502.01518, 

[41] ATLAS, ATLAS-GONF-2012-147, Search for new phenomena in monojet plus missing 
transverse momentum final states using lOfb-1 of pp collisions at sqrts=8 tev with the atlas 
detector at the Ihc, . 

[42] CMS, V. Khachatryan et. ai. Search for dark matter, extra dimensions, and unparticles in 
monojet events in proton-proton collisions at s/s = 8 TeV, arXiv: 1408.3583. 

[43] ALEPH, DELPHI, L3, OPAL, SLD, LEP Electroweak Working Group, SLD 
Electroweak Group, SLD Heavy Flavour Group, S. Schael et. ai. Precision electroweak 
measurements on the Z resonance, Phys.Rept. 427 (2006) 257-454, |hep-ex/0509008|. 

[44] Fermi-LAT, M. Ackermann et. ai. Searching for Dark Matter Annihilation from Milky Way 
Dwarf Spheroidal Galaxies with Six Years of Permi-LAT Data, arXiv: 1503.02641, 

[45] L. Goodenough and D. Hooper, Possible Evidence Eor Dark Matter Annihilation In The Inner 
Milky Way Erom The Eermi Gamma Ray Space Telescope, arXiv: 0910.2998. 

[46] D. Hooper and L. Goodenough, Dark Matter Annihilation in The Galactic Genter As Seen by 
the Eermi Gamma Ray Space Telescope, Phys.Lett. B697 (2011) 412-428, |arXiv: 1010.2752|. 

[47] D. Hooper and T. Linden, On The Origin Of The Gamma Rays Prom The Galactic Genter, 
Phys.Rev. D84 (2011) 123005, |arXiv: 1110.0006|. 

[48] K. N. Abazajian and M. Kaplinghat, Detection of a Gamma-Ray Source in the Galactic 
Genter Gonsistent with Extended Emission from Dark Matter Annihilation and Goncentrated 
Astrophysical Emission, Phys.Rev. D86 (2012) 083511, |arXiv: 1207.6047|. 

[49] T. Daylan, D. P. Finkbeiner, D. Hooper, T. Linden, S. K. N. Portillo, et. al. The 
Gharacterization of the Gamma-Ray Signal from the Gentral Milky Way: A Gompelling Gase 
for Annihilating Dark Matter, arXiv: 1402.6703, 


19 



[50] F. Galore, 1. Cholis, and C. Weniger, Background model systematics for the Fermi GeV excess, 

arXiv:1409.0042, 

[51] F. Galore, I. Gholis, G. McGabe, and G. Weniger, A Tale of Tails: Dark Matter Interpretations 
of the Fermi GeV Excess in Light of Background Model Systematics, arXiv: 1411.4647, 

[52] http://lepsusy.web.cern.ch/lepsusy/www/sleptons_summer04/slep_final.html. 

[53] M. Gamy, A. Ibarra, and S. Vogl, Signatures of Majorana dark matter with t-channel 
mediators, arXiv: 1503.01500. 

[54] AMS Collaboration, M. Agnilar et. al.. First Result from the Alpha Magnetic Spectrometer 
on the International Space Station: Precision Measurement of the Positron Fraction in 
Primary Gosmic Rays of 0.5-350 GeV, Phys.Rev.Lett. 110 (2013), no. 14 141102. 

[55] AMS Collaboration, L. Accardo et. al. High Statistics Measurement of the Positron 
Fraction in Primary Gosmic Rays of 0.5500 GeV with the Alpha Magnetic Spectrometer on 
the International Space Station, Phys.Rev.Lett. 113 (2014) 121101. 

[56] AMS Collaboration, M. Agnilar et. al.. Electron and Positron Fluxes in Primary Gosmic 
Rays Measured with the Alpha Magnetic Spectrometer on the International Space Station, 

Phys.Rev.Lett. 113 (2014) 121102. 

[57] L. Bergstrom, T. Bringmann, I. Gholis, D. Hooper, and G. Weniger, New limits on dark 
matter annihilation from AMS cosmic ray positron data, Phys.Rev.Lett. Ill (2013) 171101, 
|arXiv: 1306.3983|. 

[58] A. Ibarra, A. S. Lamperstorfer, and J. Silk, Dark matter annihilations and decays after the 
A MS-02 positron measurements, Phys.Rev. D89 (2014), no. 6 063539, |arXiv: 1309.2570|. 

[59] M. Girelli, E. Del Nobile, and P. Panel, Tools for model-independent bounds in direct dark 
matter searches, JGAP IZIO (2013) 019, |arXiv: 1307.59551. 

[60] LUX Collaboration, D. Akerib et. al. The Large Underground Xenon (LUX) Experiment, 
Nucl.Instrum.Meth. A704 (2013) 111-126, |arXiv: 1211.3788]. 

[61] E. Del Nobile, G. B. Gelmini, P. Gondolo, and J.-H. Huh, Update on Light WIMP Limits: 
LUX, lite and Light, JCAP 1403 (2014) 014, jarXiv: 1311.4247|. 


20 




